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We present a careful analysis of the temperature dependence of the nuclear spin-lattice relaxation 
rate 1/Ti in gapped quasi-lD Heisenberg antiferromagnets. It is found that in order to estimate 
the value of the gap correctly from 1/Ti the peculiar features of the dispersion curve for the triplet 
excitations must be taken into account. The temperature dependence of 1/Ti due to two-magnon 
processes, is reported for different values of the ratio r = J±/ Jy between the superexchange constants 
in a 2-leg-ladder. As an illustrative example we compare our results to the experimental findings 
for ^^Cu l/Ti in the dimerized chains and 2-leg-ladders contained in Sri4Cu2404i. 



PACS numbers: 76.60.Es, 75.40. Gb, 74.72. Jt 

The many pecuhar aspects of quasi one-dimensional 
quantum Heisenberg antiferromagnets (IDQHAF) have 
stimulated an intense research activity during the 
last decade Moreover, the recent observation 

of superconductivity in the 2-leg-ladder compound 
(Sr,Ca)i4Cu24 04i Q and the occurrence of a phase sep- 
aration in high temperature superconductors (HTSC) in 
hole-rich and hole-depleted regions analogous to spin- 
ladders have brought to a renewed interest on 
IDQHAF. One of the relevant issues is wether the spin- 
gap observed in some of these IDQHAF is related to the 
one observed in the normal state of HTSC For these 
reasons many NMR groups working on HTSC have fo- 
cused their attention on these systems and on the deter- 
mination of the spin-gap values in pure and hole-doped 
compounds [^-p^. However, since the early measure- 
ments, a clear discrepancy between the values for the 
gap (A) estimated by means of nuclear spin-lattice re- 
laxation (l/Ti) and susceptibility (or Knight shift) mea- 
surements has emerged ||^ . In many compounds the gap 
estimated by means of 1/Ti using the activated form 
1/Ti C!C exp{—A/T) derived by Troyer et al. [T^ , turned 
out to be ~ 1.5 times larger than the one estimated by us- 
ing susceptibility or inelastic neutron scattering measure- 
ments (see Tab. 1). Many attempt models, theoretical 
p8[ or phenomenological |]l^ , have tryed to explain these 
differences, however, while they were able to describe the 
findings for some compounds they were not able to ex- 
plain the results obtained in other gapped IDQHAF. In 
fact, as can be observed in Tab. 1, while for certain 2- 
leg-ladders ||l^] an agreement between the gap estimated 
from Ti and through other techniques is found, in sev- 
eral other systems it is not [^0-|2). It is interesting to 
observe that the IDQHAF where the agreement is ob- 
served are the ones in the strong coupling limit, namely 
either dimerized chains or 2-leg-ladders with a superex- 
change coupling along the rungs much larger than the one 
along the chains. Therefore, one can conclude that the 
disagreement is not always present and has to be associ- 



ated with the peculiar properties of the spin excitations in 
each system, i.e. with the form of the dispersion curve for 
the triplet excitations. In this manuscript we will show 
that the discrepancy relies essentially on the use for 1/Ti 
of an expression which is valid in general only at very low 
temperatures (T < 0.2A) and its application to higher 
temperatures depends on the form of the dispersion curve 
for the triplet spin excitations. In particular, for dimer- 
ized chains the validity of a simple activated expression 
extends to higher temperatures than for a 2-leg-ladder. 
As an illustrative example we will analyse the tempera- 
ture dependence of 1/Ti for the ^^Cu nuclei in the dimer 
chains (Cu(l)) and in the 2-leg-ladders (Cu(2)) contained 
in Sri4Cu2404i @. 

In the following we will consider the contribution to nu- 
clear relaxation arising from 2-magnon Raman processes 
only. Namely, we will assume that although the system 
is not in the very low temperature limit {T <ti A), the 
temperature is low enough (T ^ A) so that 3-magnon 
processes as well as the spin damping can be neglected. 
If the large value of the gap derived by means of l/Ti 
was due to these contributions, which are proportional to 
exp{—2A/T), one should observe some discrepancy also 
for the IDQHAF in the strong coupling limit, at vari- 
ance with the experimental findings (see Tab. 1). The 
approach we use follows exactly the same steps outlined 



in the paper by Troyer et al. |17| where, by assuming a 
quadratic dispersion for the triplet excitations (valid for 
T < A), namely 



E{k.^) = 1 + a(fc - 
in units of A, they found that 



(1) 



l/T^ ^ ^2!^_l_e-o/2Tg-A/T(o 80908 _ in{u;jT)) 



(2) 



with ujo the resonance frequency and Ao the hyperfine 
coupling constant. We remark that there is a factor 4 
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difference with respect to the equation reported by Troyer 
et al. Il^, which is related to a different definition of 
the hyperfine hamiltonian and of the dispersion curve. 
The values of the hyperfine constants are Ao — 120 kOe 
for 63Cu(2) and Ao = 29 kOe for ^^Cu{l) In 
the case of a general form for the dispersion relation, by 
considering that the low-energy processes are the ones 
corresponding to an exchanged momentum g ~ and 
q ~ ~2kx, one can write the contribution related to 2- 
magnon Raman processes in the form p7|] 



1/Ti 



7r2 knA 



dkx — ^=^= 



-E{K)/T 



2UJ, 



(3) 



where E{kx) is the dispersion relation for the triplet spin 
excitations, normalized to the gap value, whereas v{kx) = 
dE{kx)/dkx. For a 2-leg-ladder a general form describing 
E{kx) is 

E{k,)^ = E{kx = Qfcos\^) + ,sz?i2(^) + c„W(fc,) 

(4) 

which is strongly dependent on the ratio r = Jj_/J|| be- 
tween the superexchange coupling along the rungs and 
along the legs. We have taken the dispersion curves de- 



rived by Oitmaa et al. |20| from an extensive series studies 
and estimated the parameters E{kx = 0) and Co accord- 
ingly. Then, starting from Eqs. 3 and 4, by means of 
a numerical integration one can derive directly 1/Ti for 
a 2-leg-ladder for different values of r. It should be re- 
marked that for r of the order of unity the dispersion 
curve for the triplet excitations has a maximum around 
a wave-vector km (see Fig. 1) and also low-energy 
processes from fc„i — k^ to km + k^ could contribute to 
the relaxation. However, this processes should become 
relevant only at T > A, where also 3-magnon processes 
and the damping of the spin excitations become relevant. 

In Fig. 2 we report the results obtained on the basis 
of Eqs. 3 and 4 for ^^Cu(2) for different values of the 
superexchange anisotropy r. One observes that while for 
the dimerized chains, corresponding to the limit r 3> 1, 
1/Ti follows an activated behavior as the one given in 
Eq. 2, for the 2-leg-ladders with r ~ 1 one observes some 
differences with respect to the simple activated behavior 
already at temperatures T < A/4. This analysis points 
out that for a 2-leg-ladder with r of the order of unity 
it is not correct to estimate the gap from l/Ti by using 
Eq. 2, at least for T > A/4. In fact, it is noticed that 
the quadratic approximation for the dispersion curve be- 
comes valid for a more restricted range of kx around 7r/a 
as r decreases (see Fig. 1). This seems to contradict the 
results reported in Fig. 2a, where the departure from 
the quadratic approximation is found more pronounced 
for r = 1 than for r = 0.5. However, this artifact is 
related to the choice of the horizontal scale, namely to 



have reported 1/Ti vs. A/T, since A increases with r 
§. In fact, if we report l/Ti vs J\\/T (Fig. 2b), with Jy 
independent of r, one immediately notices that the devi- 
ation from the quadratic approximation starts at lower 
temperatures for the lowest value of r. 

One can then analyse the experimental data on the ba- 
sis of Eq. 3 by taking the value for the gap estimated by 
other techniques and check if there is an agreement. We 
have fit the experimental data for ^'^Cu(2) (Fig. 3b) and 
63Cu(l) (Fig.3a) in Sri4Cu2404i by taking A = 450 K 
and A = 120 K, respectively, as estimated from suscep- 
tibility or NMR shift data |7pC[| . In both cases we find a 
good agreement between theory and experiment by tak- 
ing 1 > r > 0.5 for the ladder site and r ^ 1 for the 
chain site. If the data for ^■^Cu(2) were fitted according 
to Eq. 2 one would derive a value for the gap around 650 
K, a factor 1.5 larger than the actual value (see Tab. 1). 

For r = 1 also a quantitative agreement with the ex- 
perimental data for ^^Cu(2) is found. However, this fact 
seems to be at variance with the estimates by Johnston 
PH based on the analysis of DC susceptibility data and 
with the recent findings by Imai et al. ||l^ based on the 
study of ^^O NMR shift anisotropy, where a value for 
r ~ 0.5 was derived. If we take this value for r we 
find that the experimental data are a factor ~ 8 larger 
than expected. This disagreement could originate, at 
least partially, from having considered for the = 2kx 
processes the values for the | < —kx\Sz\kx > P matrix 
elements estimated by Troyer for the case r = 1 [ p7[ . 
One has also to mention that the estimate of the hyper- 
fine coupling constants could suffer from some uncertain- 
ties, particularly the contribution from the transferred 
hyperfine interaction with the neighbouring Cu^^ spins. 
This contribution should be particularly relevant for the 
^^Cu(l) nuclei while it should be small for ^^Cu(2). How- 
ever, it must be recalled that since 1/Ti depends quadrat- 
ically on the hyperfine coupling constant even for ^^Cu(2) 
sizeable corrections can be exepected. Finally it has to be 
observed that in these systems the low-frequency diver- 
gence of 1/Ti is cut because of the finite coupling among 
the ladders (or chains), introducing another correction to 
the absolute value of 1/Ti. 

The low- frequency divergence of 1 /Ti was found to fol- 
low the logarithmic behavior reported by Troyer et al. 
|]l7| (see also Eq. 2) and does not change upon varying 
the anisotropy factor r, for r > 0. In fact, the form of 
this divergence is related to the shape of the dispersion 
curve close to kx — tt/o, where it is always correctly ap- 
proximated by a quadratic form for r > 0. 

In conclusion we have presented a careful analysis of 
the problem of estimating the spin-gap from nuclear spin- 
lattice relaxation measurements in IDQHAF. It is found 
that in order to estimate correctly the gap one should ei- 
ther perform the experiments at temperatures T ^ 0.2A 
where in many cases other contributions to the relax- 
ation process emerge ||^,^,^, or use an appropriate ex- 
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prcssion for 1/Ti which takes into account the form of 
the dispersion curve for the triplet excitations. Then a 
good agreement for the gap value estimated by means of 
1/Ti and other techniques is found, allowing also to de- 
rive information on the anisotropy of the superexchange 
constants. 

We would like to thank D. C. Johnston for useful dis- 
cussions. The research was carried out with the financial 
support of INFM and of INFN. 



FIG. 1. The dispersion curves for the triplet excitations in 
a 2-leg-ladder are reported for different values of the superex- 
change anisotropy r = Jx/J||. The dotted line shows the 
quadratic approximation of the dispersion curve for r >> 1 
(see Eq. 1). 

FIG. 2. a) Nuclear spin-lattice relaxation rate as a function 
of A/T for three different values of the ratio r = Jx/Jy. The 
solid lines give the results obtained from Eqs. 3 and 4 after a 
numerical integration for u!o/'2v = 15 MHz, while the dotted 
lines show the corresponding behavior by using the quadratic 
approximation (see Eq. 2). b) The same data as in a) are 
now reported as a function of Jy /T. 
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FIG. 3. a)Temperature dependence of '*^Cu(l) l/Ti in a 
Sri4Cu2404i single crystal for H \\ b. The solid line shows the 
behavior expected from Eq. 3 by taking a gap A = 120 K. 
b)Temperature dependence of ''■^Cu(2) l/Ti in Sri4Cu2404i 
for H \\ b. The data were obtained either in oriented powders 
(circles) or single crystals (squares). The solid line shows the 
behavior according to Eq. 3 for r = 0.5 by using the same 
value for the gap derived from '^^Cu(2) NMR shift (A = 450 
K), while the dotted line gives the corresponding behavior 
obtained using the quadratic approximation (Eq. 2) with 
A = 650 K. 



TABLE I. Values for the gap A between singlet and triplet 
excitations for different IDQHAF, estimated from 1/Ti using 
Eq. 2 and from DC susceptibility (or NMR shift) measure- 
ments. In the last column the ratio for the values of the gap 
estimated by the two techniques is reported. 
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Fig. 1 (Meizi and Carretta) 
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650 K 
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(2-leg-ladder) 
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(2-leg-ladder) 
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650 K 


660 K 
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(dimers) 
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1.25 


(S=l chain) 
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(S=l chain) 
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